Microscopic and Macroscopic Behaviors of Palatini Modified Gravity Theories 
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We show that, within modified gravity, the non- linear nature of the field equations implies that the 
usual naive averaging procedure (replacing the microscopic energy-momentum by its cosmological 
average) is invalid. We discuss then how the averaging should be performed correctly and show 
that, as a consequence, at classical level the physical masses and geodesies of particles, cosmology 
and astrophysics in Palatini modified gravity theories are all indistinguishable from the results of 
general relativity plus a cosmological constant. Palatini gravity is however a different theory from 
general relativity and predicts different internal structures of particles from the latter. On the other 
hand, and in contrast to classical particles, the electromagnetic field permeates in the space, hence 
a different averaging procedure should be applied here. We show that in general Palatini gravity 
theories would then affect the propagation of photons, thus changing the behaviour of a Universe 
dominated by radiation. Finally, Palatini theories also predict alterations to particle physics laws. 
For example, it can lead to sensitive corrections to the hydrogen energy levels, the measurements 
of which could be used to place very strong constraints on the properties of viable Palatini gravity 
theories. 

PACS numbers: 04.50.+h 
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I. INTRODUCTION 

Extensions of General Relativity (GR) have always re- 
ceived a great deal of attention. Such theories are moti- 
vated by quantum gravity models and by the wish to find 
phenomenological alternatives to the standard paradigm 
of dark matter and dark energy [ 2, [ 3 |,[j, l a , | 6 |,l7|,[8|i l 9 l . llO . 
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Alternatives to Einstein gravity generally result in non- 
linear corrections to the field equations. The application 
of these equations to macroscopic {e.g., cosmological) 
scales involves an implicit coarse-graining over the mi- 
croscopic structure of matter particles. However, when 
there are extra non- linear terms in the field equations, 
a priori, the validity of the usual coarse- gra ining proce- 
dure can no longer be taken for granted |22|. Hence, as 
discussed in ^27i . i28i] , it may be important to take into ac- 
count the microscopic structure of matter when applying 
the field equations to macroscopic scales. 

Unfortunately, up until now that has not been the 
common practice [H, [s^, lU, [s^, [111. This is proba- 
bly because, in GR as in Newtonian gravity, the micro- 
scopic structure of matter is not particularly important 
on macroscopic scales. It is standard practice to replace 
the metric, gab, and the energy momentum tensor. Tab, 
with some averages of them that coarse-grain over the 
microscopic structure of matter particles. This simple 
procedure only works, however, because on microscopic 
scales the equations of GR are approximately linear. 
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In this article we show that such an approach cannot 
simply be applied to modified gravity theories without 
a detailed analysis of the energy-momentum microstruc- 
ture. Indeed, the naively averaging over the microscopic 
structure will generally lead one to make incorrect pre- 
dictions, and inaccurate conclusions as to the validity 
of the theory [l^. Indeed, it's possible for a theory 
that deviates significantly from GR at the level of the 
microscopic field equations to be indistinguishable from 
GR when correctly coarse-grained over macroscopic {e.g., 
cosmological) scales. We illustrate this point for a class 
of modified gravity theories in which the Ricci scalar R in 
the Einstein-Hilbert action is replaced by some function 

f{R, R'^^'Rab)- 

It is well known that whenever such a replacement is 
made, the field equations for the action can be derived ac- 
cording to two inequivalent variational approaches: met- 
ric and Palatini. In the former case, Rab and R are taken 
to be constructed from the metric (jab which governs the 
conservation of the energy momentum tensor, and the 
field equations are found by minimizing the action with 
respect to variations in gab. In the alternative. Palatini 
approach, R = Rabg""^ where Rab is a function of some 
connection field F^^ that is, a priori, treated as being 
independent of gab- The field equations are then found 
by minimizing the action with respect to both F^^ and 
gab- If f{R, R'^'^Rab) = R-2A {i.e. GR with a cosmolog- 
ical constant) then the two approaches result in the same 
field equations. Otherwise they are generally different. 

Before going into further details there is one point to 
be noticed. It has been argued that Palatini approach, as 
outlined above, swaps one theory for another when ap- 
plied to f{R) actions ^43,]. Whether or not this is indeed 
the case, the Palatini f{R) field equations are mathemat- 
ically equivalent to a w = — | Brans-Dicke theory with a 
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potential, and so certainly do correspond to a mathemat- 
ically valid, and widely studied, modified gravity theory, 
even if it is not technically derivable from an /(i?) action. 

This work is organized as follows. In § |TT] we briefly 
introduce the main ingredients of Palatini modified grav- 
ity theories, derive the gravitational field equations for 
a general f {R, R"''' Rab) Lagrangian and discuss their be- 
haviour in vacuum. In § IIIII we explain why the popular 
naive averaging procedure (in which one simply replaces 
the quantities in the field equations with some coarse- 
grained averages) fails in some cases, and detail our new 
averaging method. We then reconsider the particle kine- 
matics, cosmology, astrophysics and atomic physics thor- 
oughly using the new approach and compare our results 
with old results in the literature in the subsequent sec- 
tions. In § IIVI the motion of classical particles (clumps 
of energy density in tiny patches in between which there 
is vacuum) is considered, and we find that particles in 
Palatini theories move in exactly the same ways as they 
do in GR, their active gravitational, passive gravitational 
and inertial masses are all equal, and most importantly 
the predictions on cosmology and astrophysics are also 
the same as those of GR. §|V]is devoted to an analysis of 
the behaviour of electromagnetic field in Palatini theo- 
ries: in contrast to classical particles the electromagnetic 
field permeates in the space and its averaging is a bit dif- 
ferent. We find that in general Palatini theories (albeit 
not ones where f{R,R°'^Rab) = f{R)) the propagation 
of photons is altered as compared with GR, and the uni- 
verse dominated by radiation will also behave differently. 
§ IVII then considers the atomic physics. We argue that 
for atomic physics calculations it is more convenient to 
work in the Einstein frame metric and show that the mat- 
ter Lagrangian is modified at the field theoretic level. In 
particular, the atomic energy levels now depend on the 
modification very sensitively and experimental data puts 
very strong constraints on any Palatini-type deviations 
from GR. Although the analysis is performed in the Ein- 
stein frame, we show that the resulting experimental con- 
straints are independent of one's frame choice. We finally 
summarize in § IVIIl 



noted by F^^. In what follows we shall use several differ- 
ent notations and for clarity we define them here. We use 
gab to denote the metric whose Levi-Civita connection is 
F^f, and as such we have Rab — ^ab(F) = Rabig), the 
Ricci scalar calculated from this metric is i? = g'^^Rab] 
in a similar way for the matter metric gab we have 
Rab = Rab{^) = R{g) and R = g°-^Rab- Besides this, 
we also need the mixed contractions TZ = g°'^Rab^ T^'^'' = 
gadgbcj^^^^ -^b^ ^ gbcj^^^ further define the covari- 

ant derivatives Vc and Vc to be compatible with the con- 
nections F and F respectively, i.e., "S/ cgab — cgab = 0. 
To be consistent with these conventions we shall rename 
the f{R,R''''Rab) theories as f {71,71'''' Rab) theories from 
now on, and these new notations clearly show that the 
theories at hand are neither metric ones nor pure affine 
ones. It will become clear below how these different quan- 
tities relate with each other. Note also that we assume 
Rab to be a symmetric tensor (if it contains antisymmet- 
ric part then the field equation will be spoiled). 



A. The Action 

For Palatini f {71,71°''' Rab) gravity we start from the 
following action 

Sf(n,n^^R^,) = J d^^V^j^f {71,71"" Rab) 

+ Snia.Ucr{gfj.,y,'4>i), (1) 

where n = SnG with G being the gravitational constant; 
c = h = 1. S'niattcr IS the matter action depending only 
on the matter metric gab and specific matter species ipi, 
and not on F^^. This means that the energy momentum 
tensor, defined as: 



2 SSn 



-9 sr" 



(2) 



is conserved with respect to (jab and the particle geodesies 
are determined by the metric gab- The conservation law 
is 



II. PALATINI f{R, RabR"'') THEORIES 

In this section we briefly summarize the main ingredi- 
ents of Palatini modified gravity theories. 

In general, to modify gravity one could add functions 
of the curvature invariants R, R^^Rab, R"'"'^'^ Rabcd to the 
standard Einstein-Hilbert action. In the Palatini varia- 
tional approach, the case of R""'^'^ Rabcd has not yet been 
explored up to date, so in this paper we shall focus on 
the special class of theories where / = f {R, R"" Rab) , 
the cosmology of which has recently been the subject 
of much interest. We stress that, as mentioned in § U 
in these theories the Ricci tensor Rab is constructed from 
the connection F^j, which is generally not the Levi-Civita 
connection of the matter metric gab, which is instead de- 



^"Tab = 0, 

where Voffbc = as defined above. 

B. Field Equations 

For the sake of convenience and clearness, we shall de- 
fine $ = 7^, x = 7Z''''Rab and = 7Z\ for a general 
f {71, 71"'' Rab) = f{^,x) theory. Minimizing the action 
with respect to the variations in the connection, F^^, gives 



Vx^'c) , (3) 



where /,$ = a/($,x)/5$, /,x = 5/($,x)/5x, gab is 
the metric whose Levi-Civita connection is F^^. In f{7Z) 
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theories where f,^ = 0, or in the cases when K^c oc 6^c, 
the two metrics are related conformally: gab — f.'S>9ab- 
More generally, the relationship between the two metrics 
is a dis formal one. 

We define the matrix K by = K°'i, and minimize 
the action with respect to variations in gab to find 



/,$K + 2/^K2-i/I = KT, 



(4) 



where T"b = g'^'^Tbc and I is the (4 x 4) unit matrix. $ 
and X are then given by the following algebraic relations 



$ = tr K, X = 
The trace of Eq. (j4|) reads 

/^$<i> + 2/,xX - 2/ 
where T = T° . Defining 



tr 



(5) 



(6) 



— f% 



det I 



2lx 



K 



-,1/2 



it is straightforward to check that 



and using Eq. i 

— g'^'^f.i'Rcb + '^f,x9'^'^9'^'^RebRcd 

and so Eq. ^ is equivalent to 

G\ig) = R\ig) - \Rig)S\ 

where the potential, V{'^, x)^ is given by 

/,$$ + 2/,xX-/('i>,x) 



(7) 



^('i>,x) = 



2k 



(8) 



When written in terms of gab, the Palatini field equa- 
tions, Eq. (O, are essentially those of General Relativity 
but with a modified source term given in terms of the nat- 
ural energy momentum tensor, T°'b, together with a po- 
tential term Eq. ([5]) . This correspondence between Pala- 
tini theories and modified source theories is well known. 

We could also use Eqs. ([3l HI [6|) to rewrite the Palatini 
field equations entirely in terms of what we will refer to as 
the natural, or matter, metric gab (since it is the metric 
that naturally appears in the matter action). In this case 
the field equations for a general Palatini f{TZ,Tl"^^Rab) 
theory are fairly unwieldy and so we only present the 
form they take in the special case when / = /(T?.)- Defin- 
ing F = f,^ we have: 

Gab = -^kT;;, - igafc f"^ - j + -^(VaVfo - gafcD)i^ 



haF\JbF-^gab^^F\J,F 



(9) 



where Gab = Rab — \Rgab is the Einstein tensor con- 
structed from (jab and □ = V^Va- Since this equation 
involves second order derivatives of F (or equivalently of 
T), it is generally difficult to solve while it is often eas- 
ier to solve the equations in terms of gab- Additionally, 
the curvature of gab is often much larger, even over small 
scales, than one would have expected based on the be- 
haviour of the metric in GR. Consider for instance the 
trace of Eq. jl]): 



iUF 3(VF)2 



F 



2F^ 



(10) 



Based on GR one might expect that if F = w 1, 
R ~ 0{—kT), however from Eq. pH]) it is clear that un- 
less F, which depends algebraically on T, is only varying 
very slowly, one may actually have R 3> —kT, and so 
that gravity, as described by curvature of gab, is actually 
much stronger than one would naturally expect. The im- 
mediate upshot of this is that it may not be appropriate 
to take (jab ~ Vab over laboratory scales as one might nor- 
mally expect to be possible. We discuss the important 
implications of this later in § IIVDI and return to field 
equations in terms of the metric gab in § IIVFI 



C. Behaviour in Vacuum 

The behaviour of these theories in a vacuum is, as we 
shall see later, of great importance. In vacuum, g°''^Tbc — 
0, or possibly g°''^KTbc = — AofS^t, Aq = const, in either 
case Eq. @ gives K = <i>oI/4, and so $ = $o, X = *o/4 
where $o is given by Eq. © 



-4Ao + 2/o. 



(11) 



where we defined /o = /($o, $o/4), /,<i>o = /,<t.($o, *o/4) 
and /^o = /,x($o,$o/4)- Denoting Qo = Q($o,$g/4) 
and Vo = V{<^o, *o/4), we have 



Qo = 



(2/o - 4Ao) 



and so from Eq. ([7]) 

kT % [r'^^Tbc = -Ao5% 



$0 

/o - 4Ao 



(12) 
(13) 



8(/o - 2Ao) 

-Aeff($o)J% 



5\ 



It is very well known that the vacuum field equations 
of f{TZ,Tl°'^Rab) Palatini theories are equivalent to those 
of General Relativity with an effective cosmological con- 
stant Aeff($o)- One should appreciate however that, de- 
pending on the specific form of /, (1) there may be more 
than one value of $o that satisfies Eq. (fTT|) and so the 
vacuum may not be unique in these theories; (2) the value 
of Aoff could be either positive or negative. 
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III. WHY AVERAGING MATTERS? 



particles we have instead $ = $s where 



It is important to stress that the equations we present 
in the above section are all microscopic field equations, 
which is to say that they are only certainly valid when one 
has taken into account all of the microscopic structures 
in the distribution of energy and momentum described 
by Tab- A priori there is no reason to expect these equa- 
tions to remain valid if Tah and T are replaced by some 
coarse-grained average. In particular for non-relativistic 
baryonic matter, there is no reason to believe that we can 
course grain over the peaks in density centered on each 
nuclei. In chameleon scalar field theories [s^, UHl where 
there is a strong coupling between the scalar field and 
matter, for example, the effective macroscopic field equa- 
tions do not necessarily look like the microscopic ones 
(see Ref. § IV). 

In standard General Relativity this problem does not 
emerge, because R depends linearly on p, so the average 
value of R could be calculated simply by replacing the 
microscopic value of p with its average. This averaging 
only works, however, because of the linear dependence of 
R on p. 

In the Palatini modified gravity theories however, as we 
have seen, Rab and TZ do not generally depend linearly on 
T (and hence p). Furthermore because f{TZ,TZ'^^Rab) is 
a nonlinear function of TZ and TZ°'''Rab, there is no reason 
to expect the averaged value of / in a region of space to 
be equal to the / of the averaged values of TZ and Tl°'''Rab', 
indeed this will generally not be the case. Even so, in al- 
most all of the literature the microscopic field equations 
are simply assumed to apply on macroscopical scales. For 
instance, when the cosmology of these theories is dis- 
cussed the microscopic field equations are solved with Tab 
replaced by its cosmological average. There is no reason, 
a priori, to expect an analysis conducted along these lines 
to be valid. In fact as we shall show below, when the av- 
eraging is performed properly, the actual coarse-grained 
behavior of the theory is very different from what one 
would find if the microscopic field equations were naively 
applied to macroscopic scales. 

Now let us discuss in more detail how this nai've aver- 
aging could lead to incorrect results. Consider a body, or 
a region of space, that, microscopically, contains N non- 
relativistic particles {e.g. nuclei) each with (microscopi- 
cally) uniform density pc and each occupying a volume 
(or with an average volume) Vp. The body as a whole 
is taken to have volume Vtot. The space in between the 
particles is assumed to be filled with some diffuse (non- 
relativistic) substance with microscopic energy density 
Ps {e.g. a diffuse electron cloud or, if ps = 0, just empty 
space) . 

For clarity we take the f{TZ) case as an example. By 
Eq. ([6]), inside the particle we have $ = $c such that 



/,4,($,)$c-2/($e) = 



-UPc 



where we have assumed Pc ~ const, and that the particles 
are non-relativistic. Similarly in the space between the 



/,$($«)$«- 2/($,) = kTs = -nps. 

It is now straightforward to work out the true volume av- 
eraged value of /(<&) and /,$(<!>). By volume averaging 
of a quantity Q{x) we mean 



{Q) = 

It is then clear that 



(P) 



and similarly 



1 



NVp 

Vtot 



NVp 

Vtot 



Vtot 



(/(<&)) = /($,) + [/($,) - /($, 



NVp 

Vtot 



(14) 



(15) 
(16) 



{F{^)) = /^^($,) + [/^^($,)_/^<,(<i>,)]-^, (17) 

Vtot 

and so on. In general the averaged value of the quantity 

Q{x) is 



Qs + {Qc - Qs 



Vtot 
{P) ~ Ps 



(18) 



Generally the microscopic field equations give Q as a 
function of T « — p, i.e., Q{x) = Q{p{x)) = Q{p). If 
one simply replaced p by (p) then one would think that 
{Q) — Qa = Q{{p))- It is obvious however that unless 
Q{p) — ai + a2p (where ai 2 do not dependent on the Tab 
components) (Q) ^ Qa. 

Because in the Palatini f{TZ) theories / is usually a 
nonlinear function of TZ, we conclude that, except for 
isolated values of TZ we have {f{TZ)) ^ f{{TZ)). As such 
the effective macroscopic {i.e., coarse-grained) behavior 
of Palatini f{TZ) theories will not, as is almost always 
assumed, be well described by the microscopic equations 
with Tab replaced by (Tab). A generalization of the above 
argument to the Palatini f{TZ'^^Rab) case is straightfor- 
ward if more complicated, due to the similar nonlinear 
and algebraic nature of the field equations. As a result 
the predictions made by applying the microscopic field 
equations in Palatini modified gravity theories to macro- 
scopic {e.g., astrophysical and cosmological) settings can- 
not generally be trusted. 

In the following sections we shall study the behaviors of 
classical particles and radiation fields in Palatini gravity 
theories when the averaging strategy Eq. (|18p is used, and 
compare with the results obtained by naive averaging. 
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IV. MOTION OF CLASSICAL PARTICLES 

In this section we consider the motion of a number 
of (classical) microscopic particles in Palatini modified 
gravity theories. By a classical particle, we simply mean 
some localized distribution of energy and momentum, 
i.e., for particle /, Tab only inside some world tube 
yV(7). We further require that outside the world tubes 
of the particles, gravity is weak, which means that grav- 
ity is, to a good approximation, Newtonian (or indeed 
post-Newtonian). For this to be the case, we must re- 
quire that the typical separations between the particles 
are large compared to their sizes. We also assume that 
these separations are large enough so that the particles 
are effectively collisionless. In between the particles we 
have vacuum, i.e.. Tat = and so particles only interact 
with each other gravitationally. In f{TZ) theories, we can 
relax this requirement to T = and so allow for electro- 
magnetic forces between the particles, however for the 
moment we shall not do this. 

We begin by considering a very simple, and rather ide- 
alized version of this set-up, where all the particles are 
assumed to be spherically symmetric, and furthermore, 
gravity is assumed to be weak not only outside of the 
particle world tubes but also inside them. This analysis 
was first presented in [s^] and is repeated here to serve 
as an illustration of the more general result we present 
later in this work. 



A. Spherically symmetric particles 



case is well known. Indeed if one works in terms of T°'i, 
rather than T^^, the equations that must be solved are 
just those of static, spherically General Relativity where 
outside the particle, 

f\ = -A,ff(<I>o)(5%. 

Taking the following ansatz for the metric functions 

= Wir)e^^'^'-\ = 1/M^(r), 

then the tt and rr components of the Einstein equations 
are 

2dC \ , 

The solutions are [H 



W{r) = ,_ ^GMir) _ AM) ^, 



C{r) = dr'/ P + P^+P^ > 



where C is a constant of integration and 

Ai^-y($,x)-Q§K)/Q'(*,x). 

Outside the particle C is a constant given by 



(20) 



r 3 

(22) 



1. A single Particle 

We begin by considering a single, isolated, static and 
spherically symmetric particle centered at r = and has 
a radius Tp. Outside the particle T°fc = and without 
loss of generality we work in coordinates where 



gabdx^da;^ = ds^ 



=^(-^*)dt2 + e^('-'*)dr2-f-rW. 



Because the particle is static, so = 0, and the spherical 
symmetry imposes the conditions e — — Pi~Pa/'2, 
say, and — Tg — 0. We define T'V = pA + pi, 
T'^o — ^P- One could find the metric gab both inside 
and outside such a particle by directly solving Eq. 
however, it is much simpler to find gat via Eq. (O, and 
then if necessary one can recover the results in terms of 
gab by some transformations, as we shall show in § IIVFI 
Note that Eq. ([7]) is nothing but the standard Einstein 
equation with a modified energy momentum tensor: 



Ta 
h 



(19) 



^-2 Jo ^''Q2(<i>,x)m-') ' 

and for convenience we can set this constant to be zero 
by redefining the time coordinate. Then 

C'^~- I dr'r' P + P^+P^ 



2 Jo ■ ■ QH^,x)W{r') 
and so finally C can be written as 

C = r dr'r' P + P^+P^ 



(23) 



Additionally, outside the particle p = AV = and so 
M{r) — Mp — M{rp) — const. The angular compo- 
nents of Eq. ([7]) provide us the analogue of the Tolman- 
Oppenheimer-Volkoff (TOV) equation: 

dPeff , P + Pi+PA^, . ^PA 

+ Y{r) = ~3— , (24) 



S2 



where 



We can see from this expression that Tq = if and only if 
Tq — 0, and that all the symmetries of T°f, are preserved 
by T'^b- The solution for the metric functions in this 



dr Q^rW{r) 



Peft^{pi-PA-AV)/Q' 



Y{r) = (47rGPeffr' + GM - AesrVS). 
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This last equation implies that PeS is Cq continuous at 
r — Tp, and so Pcsii^p) = 0- The key thing to note is 
that, no matter what the internal structure of the par- 
ticle is, and no matter what the details of the Palatini 
theory are, outside such a particle the metric is simply 
Schwarzschild-de-Sitter, just as it is in GR with cosmo- 
logical constant Acg. This is of course precisely what one 
should expect given that we have assumed spherical sym- 
metry and the vacuum field equations are the same as GR 
with A = Aoff. As expected as it might be, this simple 
observations carries with it an important corollary which 
should be stressed: to the outside world the gravitational 
field of the particle has precisely the same form as it does 
in GR. The gravitational mass, Mp, depends on the com- 
ponents of T^b in a slightly more complicated manner 
than it does in GR. However, these components cannot 
be measured gravitationally by an external observer, and 
only Mp can be measured. Indeed, as we shall show be- 
low, Mp is, in all the usual senses, the physical mass 
of the particle, i.e., it is clearly the active gravitational 
mass but it is also equal to the passive gravitational and 
inertial masses of the particle. 



2. Multiple Particles 

Above we reviewed the spacetime about a single spher- 
ically symmetric particle in Palatini theories, and noted 
that to an external observer it was, in all theories, simply 
Schwarzschild-de-Sitter, just as in GR with a cosmologi- 
cal constant. This would be true of any modified source 
gravity where a vacuum form for the natural energy mo- 
mentum tensor, T^, corresponds to a vacuum form for 
the modified energy momentum tensor, T";,, i.e., 
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We now make a fairly straightforward generalization to 
the case of multiple particles. We work in the weak field 
limit and assume that the relative motions of the particles 
are non-relativistic, which means V ^ where d^V 

can be any one of Kr%/Q^ nT^^jQ?, KAy($,x)/Q^ 
0(e), ~ o(e) where d'^V^ = nT%/Q^ - o(e) where e 
is a small parameter. 

Correspondingly we linearized the metric gab as 

gab = Vab + hab 

with /iQo, hij 0{e) and /iq; ~ The assumption 

that the system is non-relativistic is equivalent to assum- 
ing that it is quasi-static, that is, for any quantity A each 
time derivative is suppressed by a positive power of e rel- 
ative to each spatial derivative. 

The metric inside a particle (labelled K), which would, 
by assumption, be spherically symmetric in the absence 
of any other matter fields, centered at X(^x) (t) can then be 
calculated straightforwardly by linearizing the Einstein 
equation for gab, Eq. ([7]). 



We have chosen to linearize gab rather than (jab, be- 
cause gab essentially obeys the Einstein equation, though 
with a modified source, and so provided this source is 
suitably weak the deviations of gab from rjab will also be 
small; depending on the details of the theory the same 
may not be true of g^^ (see Sections lIVDl lIVFp . We 
shall neglect the effects of the effective cosmological con- 
stant, Acff ($0), since on scales that are very much smaller 
than the cosmological horizon, it is known to have negli- 
gible effect on particle motions. 

To 0(e) the "0 component of Eq. ([7]) becomes: 

1 _ n (T"o - T\ + 2AV) 

2 - 

where, inside the if'^ particle we define T'^o = —p^^^ 
and to leading order 

^ - \Pl ^2 A )°3^0^A y(K)y{K)r{K)- 

Here y^^^ = - 2:(^)(i) and t^k) = \Jv\k)V\k) 
radial coordinate centered at (i) • Defining 25 = /iqo 
we then have: 



p + ipi-2AV{<^,x)]. (25) 



Both inside the K^^^ and in the vacuum between the par- 
ticles, Eq. ([^5]) has the solution: 

GM(^K){\x~X(^K){t)\) 



S{x\t) 



\x-X^K)\ 
-Ci^K){\x - X(^K){t)\) + U(K){x\t), 



where 



X ^ GMi j\ 



(26) 



where C(i^:)(r(/f-)) and M(x')(r(^)) are given by Eqs. 
and (I22p with r r(x) , p —* p^^^ , Pi p\^^ and pA 
p^-*. Mf^j) is the total mass of the J*'' particle and is 
also given by Eq. (j22|) with the appropriate substitutions 
in the limit ?' — > 00. 

In a similar manner, from component of Eq. ([7]) we 
have, h,j = 2*(5y where to 0(e) *,fcfc = K.f°o/2: 



GM(r(K)) 
r(K) 



(27) 



where for a particle with radius Vp 

so outside the particles D = Q. Finally for joi we fix the 
gauge so that: 



4*,o, 
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and then the Einstein equations give 

hoi^kk — — 2kT*o 

We can now write down the metric outside aU but the 
K^'^ particle: 

gabdx'^dx'' 

= - l-^^^MhI^A + 2C(r(K))-2U{x\t)\dt'{28) 

L ^K) J 



, . ^ 2GM(r(K)) 



dx^ 



Because outside the particles we have C = -D = 0, so the 
metric is given simply as 

gabdx^'dx'' = - (1- 2*A,)dt2 + (1 + 2^'jv)dx^ (29) 

where 

*jv = C/(A-) + GM(K)/r(K) = GM(j)/r(jy 

J 

3. Motion of Particles 
It is clear from Eq. (O that if we define 

we have 

Vaf"'' = 0, (30) 

where V aQbc = and we used the fact that WaG'^^{g) ~ 
with Gabig) the Einstein tensor of gab- After a lengthy 
manipulation of the Einstein equations its is well known 
that one can extract the following relation [i^ 



V^Vaf"'' = da [T"'' + t"''] = 0, 

in which T"'' = y^f"'^ 



(31) 



2k 



with 



yUb 

L 



gg U ad^ be ^ ab^ cd) J 



c j^d 



and 



with 



(32) 



n 



acbd 



lb „cd „bc„ad 

g -g g , 



^ab 



nb 



t°b is called the pseudo-tensor density of gravitational en- 
ergy and momentum and it is not a true tensor density. 
In contrast T"b is a true tensor density. Eq. ([3T|) is not a 
covariant equation because t"^ is not a tensor density and 
that it involves the ordinary rather than tensorial diver- 
gence. However, the ordinary divergence will be useful 
to obtain (for finite systems) the analogues of conserva- 
tion of energy and momentum as in classical dynamics 
by integration over some large enough space region. 

Returning to Eq. ([?!]) ■ suppose that the fC'^ particle 
is in a finite region of space, Vk, (which is a reason- 
able assumption) that contains only the K^^ particle and 
which is enclosed by a boundary, E^, located in the sur- 
rounding empty space at a sufficient distance so that in 
some open set of points including the boundary gravity 
is weak. We can then choose coordinates {t, a;*) such that 
gab = Vab + hab whcrc hab is Small. With these coordi- 
nates we have from Eq. PT|) : 

^ / d^x [T°" + 1°°] = - / d^x [T'" + 1*''] . .(33) 
dt Jvk 

Defining J^"^) = ^^^^ d^x [T°° + t°°] we have 
dJ'' 



dt 



(34) 



Using Eq. ((32|) we see that J^°^^ is itself given by a surface 
integral: 



(35) 



Using Eq. (^5]) we find that outside the particles to 0{e) 
(remember that g"'' = g"^^\/—g) 

g0" = -(l + 4*), g'^^S'i, 



g" = - 


so 








M^ K ) - 


li 


M(K)- 





-JiK) 



U,d'^, 

Similarly we find after some manipulation 



K) 



dc'xdtT^ox' 



V,. 



and since Tq depends on and t only in the combination 



(t) I at leading order we have 



P — 



Xq 



d'f% 
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So defining tt^^j = whicii is the 4-velocity of 

the center of mass of the particle we have: 

^fx) = ^V)"(V)- 

Now by calculating on '^.(^fc) and using Eq. ([34| find 
the evolution of J|^^ . The component of Eq. ([M)) gives 

M(K) = 0, 

and by noting that in the vacuum regions 
the i components gives 

= - [ d^a;f%[/,, (36) 

where the last line follows from expanding U^i about — 
•^(K)' ^-kk — inside the body, and T"o being a function 
of r(K) only at leading order. Thus to leading order in e 
we have 

= UA^U)), (37) 

which is precisely the Newtonian equation of motion, and 
precisely the same result that we find in General Rela- 
tivity. Since gab is the metric that appears in the matter 
action, it is clear that small particles will move along 
geodesies in gat- Now Eq. ([57]) tells us that small par- 
ticles will also move along geodesies in gab- Since inside 
the body, in general, gab and gab are only related disfor- 
mally, and hence will have different geodesies, this result 
may seem rather counter-intuitive. It can, however, be 
understood in a fairly simple fashion. Although the two 
metrics are in general related disformally, we noted in 
Section [II CI that in vacuum the two metrics are related 
by a constant conformal factor, i.e., up to a rescaling of 
coordinates they are the same metric. Thus outside the 
particles, the geodesies of the two metrics are the same. 
Inside a particle then the only differences between the 
geodesies of the two metrics are due to essentially local 
effects, i.e., they are due only to the local matter content. 
Such local differences cannot lead to the particle devel- 
oping an overall acceleration, as this would constitute a 
self-acceleration and hence a violation of energy and mo- 
mentum conservation. Both the natural, T^h, and the 
effective, T°'b, are conserved (with respect to gab and gab 
respectively). They also both vanish outside of the parti- 
cle implying that there is no flux of energy or momentum 
in or out of the particle. Thus the total energy and mo- 
mentum inside the particle must be conserved, and there 
can be no self-accelerations. It should therefore come as 
no surprise that particles move along geodesies in gab- 



B. Generalization 

We have seen, as we first reported in |37j], that the 
motion of spherically symmetric, classical particles in a 
general Palatini f{TZ,TZ°'^Rab) theory is observationally 
indistinguishable from the motion of the same set of par- 
ticles in pure General Relativity with a cosmological con- 
stant, Aoff. The reason for this is remarkably simple, al- 
beit unappreciated up to this point, and follows from a 
couple of very well known phenomenon, both illustrated 
above. The first is that the field equations of all Palatini 
theories are equivalent to those of GR with some modi- 
fied source T"b, when this is just an effective cosmological 
constant. The second, is that, quite remarkably, GR is in 
many ways holographic and specifically the equations of 
motions for a localised distribution of energy and momen- 
tum surrounded by vacuum can be derived by consider- 
ing surface, rather than volume, integrals over curvature 
components [ssl, [s^ . Indeed in General Relativity, even 
when a modified source is present, the gravitational mass 
of an isolated particle, as well as higher mass moments, 
can all be defined in terms of surface integrals outside the 
body and hence identified with parameters in the general 
vacuum solution (39| . 

In Ref. [1^, via a computational tour-de- force, the 
Post-2-Newtonian equations of motion for a set of N clas- 
sical particles (with arbitrary internal structure and no 
assumed symmetry) where calculated using surface inte- 
grals. The motion of the particles was found to depend 
entirely on quantities defined outside the particles which 
are naturally interpreted as the different mass moments 
of the particles. Indeed the spherically symmetric par- 
ticle analysis given above follows from a spherical case 
of the calculation done in Ref. [3^ , and so the results of 
Ref. [3^ imply that all our conclusions apply equally well 
to particles with no symmetry. 

On the largest scales, the higher mass moments play a 
relativity insignificant role and the motion of a set of non- 
relativistic classical particles is determined, to a good 
approximation, entirely by their initial positions, veloc- 
ities and their gravitational masses. As should come of 
no surprise then, one does not need to know the inter- 
nal structure of a particle to know how it moves, one 
does not, for instance, need to know the precise molecu- 
lar or atomic structure of a clump of particle to predict 
its motion under gravitational and other external forces. 
This is true for any theory that is equivalent with Gen- 
eral Relativity with a modified source provided that a 
vacuum form i.e., oc 5°'b, in the natural energy momen- 
tum tensor, T°'b, corresponds to a vacuum form in the 
effective one, T°'b- In all these theories then the modifi- 
cation of the source term does not ultimately matter as 
classical particles i.e. clumps of matter surrounded by 
vacuum, still move as particles do in General Relativity. 
Particle motions are, in so far as particle motions un- 
der gravity and other external forces go, observationally 
indistinguishable. 
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C. The Physical Mass of a Classical Particle 

In this section we consider the meaning of the physical 
mass of a classical particle. In principle, there can be a 
number of different quantities associated with a particle 
that will, in different situations, play the role of its mass. 
Firstly, we have the active gravitational mass which is 
a measure of the strength of the gravitational field in- 
duced by a body. Secondly, there is the passive gravita- 
tional mass, which determines the force that a body feels 
within a given gravitational field. Lastly, there is the in- 
ertial mass of a body, which determines how quickly a 
particle's momentum changes when a force is applied. In 
General Relativity the latter two are manifestly equal. 
Furthermore, any theory in which the passive gravita- 
tional and inertial masses are equal is said to satisfy the 
weak equivalence principle, and in the absence of any 
non-gravitational external forces the trajectory of par- 
ticles will depend only on their position and velocity and 
not on their composition. Additionally in the Newtonian 
limit of General Relativity, the two types of gravitational 
mass are also equal [4l[ , and any violations of this equal- 
ity are tightly constrained by experiments [4l| . 

We saw above the Newtonian limit of the equation of 
motion for a classical, non-relativistic, particle in both 
General Relativity and Palatini theories is the same. Our 
analysis was for spherically symmetric particles although 
the results of Ref. readily extend this to all classical 
particles. 

In general if one defines the center of mass, of 
the classical particle, inside which gravity is weak, with 
effective energy momentum tensor T"^ in the standard 
way, we have: 



where 



p 



Mp = - 



d^x'^f%. 



(38) 

(39) 
(40) 



where Vp is the volume occupied by the particle. Then 

dr" 



J" = M, 



" dt ' 

where x^^ = {t, x*^). In the absence of non-gravitational 
external forces, i.e., V aT'^b = 0, we have Mp = and: 

where U^i is the Newtonian potential due to other parti- 
cles: 



(41) 



where r(/) is the distance from x\..^ to the center of mass 
of the and M^-^ is that particle's mass. It is clear 
from this relation that Mp and the are the active 

gravitational mass of their respective particles. It is also 
clear that 



2^cm 



dt^ 



(42) 



Thus the trajectory of a classical particle is independent 
of its mass and hence its composition, and hence two 
particles at the same point would feel the same accelera- 
tion. If the difference in the acceleration of two bodies at 
the same point vanishes in one frame with one choice of 
coordinates then it obviously must vanish in all frames. 
Furthermore it is independent of which of the metrics, 
gab or gab^ one is working with when performing the cal- 
culation [1]. The absence of any differential acceleration 
is precisely the statement of the weak equivalence prin- 
ciple which in turn is equivalent to saying that the pas- 
sive gravitational mass and inertial mass of any body are 
equal. 

We should not be surprised by this. The weak equiva- 
lence principle certainly holds in GR, irrespective of the 
composition or internal structure of the particles one con- 
siders, and so it must also hold in any theory which is 
equivalent to GR up to a modified source and in which 
a vacuum form for the natural energy momentum ten- 
sor, T°b, corresponds to a vacuum form for the modified 
source term, T°h. Since, in so far as particle motions go, 
the only difference between such a theory and GR are in 
the internal composition of the particles. 

If non-gravitational forces are present and associated 
with an effective energy momentum tensor Tj'', then gen- 
erally Vaf"^ = -VaTf = f''. Assuming that the gravi- 
tational field is still dominated by the matter in the par- 
ticles, we find by repeating the analysis of Section HVAI 



dr 
"dT 



Mr, 



d^x' 



dt2 



MpU.^ 



(43) 



MpU,, 



Thus is then identified as the total, non-gravitational, 
external force on the system. It is clear from Eq. ([^5]) 
that Mp is, in addition to being the active gravitational 



[1] In fact one could apply V^T^j, = to a classical particle which is 
in static configuration, and show that the TOV equation is equiv- 
alent to the statement that the particle has no self-acceleration. 
If one forgets the internal structure (pressure gradients etc.) of 
the particle, then the gradient of <I> inside the particle would not 
be balanced and would appear in the geodesic equation of the 
particle. As # depends on the local energy density of the par- 
ticle, this would lead to the conclusion that the particle feels a 
self-force that depends on its energy density or materials and so 
WEP is violated. However, with TOV equation taken into ac- 
count, one can find that the pressure gradients and ^ gradients 
cancel, leaving no self-force on the particle. 
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mass, the inertial mass of the particle. Thus we have that 
inertial and both types of gravitational mass are equal 
to Mp in Palatini theories. In all the usual senses, i.e., 
gravitational and inertial, then Mp is the physical mass 
of the particle. However Mp is not, in general, equal to 
the particle mass in the matter action, Sm{^^ ,gab). If 
it where then Alp would depend only on T'^j,, but Mp is 
given by Eq. (j22p . and this manifestly depends not only 
on T^b but also on x) and Q($, x). 



D. The Physical Energy Momentum Tensor and 
Metric 



In the study of Palatini theories, it is common prac- 
tice to refer to T°f, as the 'physical' energy momentum 
tensor, and gab as the 'physical' metric. This is purely 
a convention inherited from GR: because if one selects a 
local inertial frame defined at some point with respect to 
(jab, then in a region around this point physics is well- 
described by the special relativistic limit of the matter 
action, i.e., 5'm(^'% ??afc)- The length of this region is es- 
sentially equivalent to the length scale over which the 
spacetime appears to be flat. 

In a general scalar tensor theory, it is usually the case 
that although gravity is modified, its strength is still 
roughly the same and hence the length scale of the cur- 
vature of spacetime is no smaller than it usually is. As- 
suming that gab ~ ^ab in a laboratory is therefore no 
more or less valid that it is in General Relativity. Un- 
der these circumstances attaching the label 'physical' to 
T^^i, and (jab seems reasonable. In Palatini theories (and 
as — > —3/2 in Brans-Dicke theories), however, things 
are different. Here, if we treat tjab as the physical met- 
ric, then gravity is much stronger over very small scales. 
Indeed one could view these theories as containing an 
additional component to the gravitational force which is 
infinitely strong but has zero range. The upshot of this is 
that the curvature of gab is much larger than one might 
naively expect, and so the length scale over which one 
can treat tjab ~ Tjab is much smaller. Indeed, depending 
on the composition of the particles considered, it may 
even be smaller than the spatial extent of the particles 
themselves. 

The presence of a new strong component to the gravi- 
tational force means that even over laboratory scales we 
can not be sure that physics will be well described by 
Smi"^^ , Tjab) ■ Attaching the label 'physical' to T"-i, and 
gab is therefore misleading. Indeed we saw above that in 
the Newtonian limit, the physical mass of a particle was 
not given by a volume integral over p = —T^a but by a 
volume integral over —T°'b. 

It is, arguably, more straightforward to treat gab and 
T"f, as being 'physical'. The gravitational side of the 
theory is then simple General Relativity as evidenced by 
Eq. Q. 

In this case we could make the definition Smi^^ , gab) — 
Smi^^, (jab), where are some redefinitions of the origi- 



nal matter fields. Since now the (new) matter Lagrangian 
depends on the metric gab whose Einstein tensor behaves 
exactly like in GR [cf. Eq. ([7])] , and on laboratory scales 
we are justified in taking gab ~ ijab, so on these scales 
special relativity applies as in GR and gravity plays neg- 
ligible role in microscopic physics. However, due to the 
redefinition of matter fields, the microscopic physics {e.g., 
field theoretic) itself is now described by some modified 
action Sm, which will generally include new interactions 
between fundamental particles. 

Throughout we have endeavored not to ascribe the la- 
bel 'physical' to either metric or energy-momentum ten- 
sor, although in Palatini theories one should be aware 
that in the frame of a laboratory here on Earth it is gab 
and T°'b, rather than gab and T°-b which behave as we 
would expect given our initiation for how things work 
in General Relativity. Ultimately, all truly measurable, 
and hence physical, quantities should be independent of 
which names one gives to which metrics or which frame 
one works in, and so the names which one gives to the 
metrics and energy momentum tensors should only be 
seen as a guide to intuition rather than having any deeper 
meaning. 



E. Coarse-graining the Energy Momentum tensor 
of Particles 



We now consider the coarse-grained form of the effec- 
tive energy momentum tensor, T'^b when microscopically 
matter is clumped into classical particles surrounded by 
a vacuum. Given that we have seen that particle motions 
in GR and Palatini theories are the same, we should ex- 
pect the coarse-grained effective energy momentum to 
have the form of the energy momentum tensor for coUi- 
sionless dust, as it would in General Relativity. We show 
this explicitly below. It should be stressed that the sim- 
ple calculation present below is not new in the context of 
general relativity [lOl , however its consequences for mod- 
ified source theories have not been appreciated so far. 

Defining t'^^ = yf^f"^ + Aesg"'' we have 

We now average T"'' over a region with fixed volume V 
and surface S. E is chosen to lie in the vacuum region 
between the particles such that on E we have T°fc — 0. 
After some algebra we find 



1 d^ 
2dt2 
1 



d^xx^x^T' 



00 



(45) 



-/d^x 



a;^x-''(r2rft^'*),o + 2x(T£f^ 



For non-relativistic particles, in a weak gravitational field 
we have - ©(« 
a factor of e^^^ so 



we have Fj^^ ~ 0(e) and each time derivative introduces 



-gr, = -AeffVJ', + 0(e) / d^xT' 



hOO 



(46) 
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Indeed using the metric found in Section [IVI and assum- 
ing that the internal structure of the particles is in equi- 
librium, Eq. ((45|) gives to 0(e) for particles with masses 



M(^) and centers of mass x^^^ (t) : 



(K) 



(47) 



Where the sum is over all of the particles inside the vol- 
ume V. Similarly, one can show that 



(48) 



[K) 



This is precisely the same as what one finds in general 
relativity, which is unsurprising given that, up to a mod- 
ified source, the two theories are equivalent, and we have 
shown that the modified source does not affect the mo- 
tion of classical particles. It is straightforward to show 
that the same in true in a cosmological background pro- 
vided the relative peculiar velocities of the particles are 
non-relativistic and the peculiar gravitational potential is 
also small. Thus it is clear that the coarse-grained energy 
momentum tensor, T^fc, is that of coUisionless dust with 
a cosmological constant Acff . On cosmological scales the 
former is known to be well approximated by a pressure- 
less perfect fluid energy momentum tensor. 



F. Changing to the Matter Metric 

In this subsection we consider the f{TV) gravity theories 
using the 'natural' metric gab and show that the results 
agree with those obtained above using the metric gab and 
with those in the literature. The method used here is 
however different from the literature - we start from the 
solutions already obtained in § IIV A II and make some 
coordinate transformation to arrive at our new results. 
This analysis is relevant because one may want to find 
the behavior of the theory in terms of the original matter 
metric. A similar analysis could be carried out for the 
general f{TZ,TZ°'''Rab) theory but that is too complicated 
and beyond the scope of this work. 

For simplicity and ease of comparison with the results 
in the literature let us for now define 



exp(^) = M^exp(2C), exp(B) 



1 

W 



(49) 



and again write the metric of the spacetime inside and 
outside a static and spherically symmetric particle as 

gabdx'^dx'' ^ -e^^'^Ut^ + e^'^'-'>dr^ +r^dn^ (50) 

where A, B are functions of the radial coordinate r due 
to the symmetry and are given through Eqs. ((49l [20]) . 

We already know that the matter metric gab is related 
with the metric gab through gab = Fgab so that the line 



element is given as 

ds^ = Q„hdx°'dx^ 



gabO 
1 



^dr^ + r^dn^] 



(51) 



On the other hand, if we start from the metric gab at 
the very beginning, then the spherical symmetry simply 
requires that the metric takes a similar form as Eq. (|50|). 
that is 



dr 



(52) 



where now the coordinates are all barred as they may be 
different from the unbarred ones in Eq. (j50p . 

It is not difficult to find out the relations between these 
two sets of coordinates, the main observation being that 
Eqs. ([51] [52]) should be equivalent to each other. Because 
there is no dependence on the angular coordinates, we 
could set, from the comparison between Eqs. ([STl [5^ . 
dQ'^ = dCl'^ and thus 



r 



Similarly we set dt — dt. Then we just need to use 



F 



A(f) 



F 



(53) 



(54) 



and Eq. ((53|l to find out the relations between A, A and 
B, B and in this way calculate explicitly the matter met- 
ric in Eq. (|52| . 

To do this, note that F is a function of T and thus of 
r. Equivalently it could also be expressed as a function 
of f. This is because there is a relationship between r 
and f by Eq. ([55]) . i.e., we could write f(r) or r(f). So 
once the form of f{TV} is known and F{r) solved as in 
Sec. IIV A 1[ Eq. ([55)1 can be used to find out r(f) and 
f(r), then F{f) — F[r{f)] could be calculated, at least 
numerically. 

In what follows we shall use a prime (star) to denote 
the derivative with respect to r (f), i.e 
F* = Then Eq. ([53]) can be written as 



^ and 



r = y/F{F)f ^ dr ^ ^jF{f){l -f -i)df 
where 7 = rF* /2F, and from Eq. ((54|) we finally obtain 
1 



exp [A(r)] 



F{f) 



exp 



exp [B{f)\ = exp B ( ^JF{f) 



A[y/F(f)7 

(1 + 7)'- (55) 



There are several points to be noted about these results: 

1. We manage to calculate the metric Eq. ((52)) with- 
out explicitly solving the complicated modified Ein- 
stein equation Eq. ([9]) . This method, when general- 
ized appropriately, should be very useful when one 
deals with the f{TZ°'^Rab) gravity theories, in which 
case, as we discussed in Sec. HIl the modified Ein- 
stein equation calculated with the physical metric 
(Jab should be very complicated. 
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2. Because F* is involved in Eq. ((55|) . i? and thus 
the metric could be discontinuous even though 
5^1/ is continuous. This discontinuity happens when 
there is a sudden change of energy density distri- 
bution, for example from p = pa — const, inside 
the particle to p = outside it. Other cases when 
there will be singularity can be found in [3ll| . Fur- 
thermore, if the size of a particle with some fixed 
mass is tiny, then 7 could be very large, making the 
metric gab deviate from 77^6 significantly. 

3. Outside the particle where F = const., for some 
choices of /(7^), e.g. f{n) = 7^" with a > 1, 7^ = 
is a vacuum solution to Eq. ^ and so outside 
the particle = 1,7 = 0. In this case obviously 
f = r and A = A^B = B outside the particle so 
that the spacetime there is exactly Schwarzschild. 
In other cases 7^ 1,7 7^ 0, the outside space- 
time is Schwarzschild-de-Sitter. Consequently the 
conclusion that the spacetime outside a particle is 
Schwarzschild-de-Sitter also holds after changing to 
the matter metric cjab- Indeed the claim in Sec. lIV Cl 
that for f{TZ,TZ°'''Rab) gravity theories the active 
gravitational mass of a particle is equal to its in- 
ertial mass is correct no matter we use gab or gab, 
which is as expected because these two are ulti- 
mately equivalent. 

As a further example, we could use the above method 
to obtain explicitly the equations that A{f), B{f) must 
satisfy. From the second equation in Eq. (I55p we have 

B{f) = B(r) + 21og(l + 7) 

the derivative of which with respect to f gives 

2 . 



B* = B'^ 



dr 1 -|- 7 



7 



dr 



F' 
F 



1 + 7 



1 



1 + 7 
1 



e'' . T.N 1 - 
— {Kp+V)r + ^^ 

b r 



1 + 7 



rF* 



p* 
) 

' F 



3f fF* 



(56) 



where in the second step we have used the relations 



7 = fF*/2F, dr = y/F{f){l + -f)df 



and 



dr 

and in the third step we have used = ^ and Eq. (PO)) . 
Similarly from the first of Eq. (j55p we have 

A(f) = A{r) - \ogF 



the derivative of which with respect to r gives 

A* A'dr F* 

A* = A' 

dr F 



1 



1 + 7 
1 



' - 1 

— h -prii^Pi + i^PA - V)r 

r r 



1 + 7 



3f fF* 
T \~F 



F 



(57) 



in which we have used 



A' 



dr 



W ' W 

and Eq. (|20p . The coupled differential equations Eqs. 
[57)1 govern the f dependence of ^, B, and by solving them 
with appropriate boundary conditions we can get A^ B. 
This is just another method to solve for A,B. It is not 
difhcult to verify that these equations are equivalent to 
Eqs. (21, 22) listed in [3^ which are derived using the 
full modified Einstein equation with the matter metric 
(jab, Eq. ([5]). This again shows the equivalence between 
these two methods. A similar set of differential equations 
could be obtained for the case of Palatini fiW^^Rab) as 
well but that is beyond the scope of this work. 

G. Discussion and Summary 

Having directly considered both the motion of classical 
particles and the coarse-graining of the energy momen- 
tum tensor for matter consisting of classical particles, we 
are now ready to discuss how the coarse-grained aver- 
aging in the Palatini modified gravity theories leads to 
effects which are distinct from what have been claimed 
using the nai've averaging. 

On the microscopic level, the matter in our Universe, 
whether it is dark or baryonic (radiation will be discussed 
in the next section), is made up of small particles. Con- 
sequently our previous analysis of the motion of particles 
in Palatini modified gravity theories is directly applicable 
to this setting. As we have seen, the effect of the Palatini 
modification to General Relativity is a change of the in- 
ternal configuration and structure of the particle, while 
outside the particle the spacetime is Schwarzschild-de- 
Sitter with a cosmological constant Aoff determined by 
the model itself. The motion of these particles in Pala- 
tini theories follows the geodesies which is precisely the 
same as in general relativity with A = Aoff. There are 
no new extra dynamical degrees of freedom in the Pala- 
tini theories, which in turn means that there is no new 
long-range forces. Any new effective force must be non- 
dynamical and act only at points; so it is entirely local 
and cannot be felt inter-particles. The closest analogue 
to this in particle physics would be Fermi's original pro- 
posal for a theory of the weak force. 

We could also see this from the viewpoint of averag- 
ing. We showed above in Section IIVEI that if T°-b de- 
scribes matter clumped into classical particles in some 
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vacuum then, when properly coarse-grained over scales 
much larger than the inter-particle separation, the mod- 
ified source term, T'^t, has the form of the energy mo- 
mentum tensor for coUisionless dust plus some effective 
cosmological constant Aoff. This is entirely as one would 
expect given our analysis of particle motions in Palatini 
theories. For N particles each with physical mass mp in 
a volume Vtot , the effective energy density of the dust is 

Pmattcr = NrUp /Vtot- 

Now we compare our above conclusion with that one 
would make when using the naive averaging by consider- 
ing some specific situations (for simplicity we take f(TZ) 
theory as an example, the case of f{Tl"'''Rab) is similar). 

1. Consider the cosmological setting. In the literature 
it has been extensively claimed that if f{TV) is cho- 
sen so that the deviation from f{Tl) = TZ grows at 
small values of TZ, e.g., f{TZ) = TZ — then the 
model could lead to a phase of accelerating expan- 
sion of the Universe at late times which is differ- 
ent from that predicted by general relativity plus 
a cosmological constant. There is nothing wrong 
with this if we can model the matter distribution 
as smooth even at the smallest scales. Since in this 
case both F(<i>) and V{^) vary as <& varies with 
time so that from Eq. ([T]) we could expect the evo- 
lution to be different from ACDM. However, as we 
have emphasized several times, it is more realistic 
to model the matter as made up of small particles 
on microscopic scales and only becoming a fluid on 
large scales after coarse-grained averaging. What 
happens after this averaging? Take F for example, 
according to Eq. we have 

{F) = Fo + [F(p,)-Fo]^ 

not 

where Fq = F{pq) = F{0). As Vp < Vtot it can be 
shown that [F{pc) — Fq] <C Fq because F does 
not depend on p linearly [F{pc) ^ Fq ^ so 
essentially we have (F) ~ Fq. Similarly {V) = Vq. 
Thus according to Eq. ([7]) the model is indistin- 
guishable from ACDM. 

2. In astrophysical environments such as the Solar 
System, the matter density is so low that NVp <^ 
Vtot- Again our analysis for the cosmological set- 
ting applies here, namely the Palatini theories 
behave indistinguishably from General Relativity 
plus a cosmological constant after coarse-grained 
averaging. This in particular means that the 
Parametrized Post-Newtonian (PPN) parameters 
we measure should be the same as those in Gen- 
eral Relativity. Note however that in some excep- 
tionally high-density astrophysical systems, such as 
neutron stars, we have NVp ^ Vtot and so the naive 
averaging may give a reasonably good description; 
in this case one could expect the model predictions 
of Palatini theories and general relativity to be dif- 
ferent. 



3. As a last example, consider the matter metric gat 
in Eq. ([5^ . Eq. ([55]) tells us that the metric func- 
tion A depends on the local value of F. Now sup- 
pose we use the naive averaging. Since the function 
A determines the results of the Rebka-Pound ex- 
periment, we can do two such experiments, one in 
the normal atmosphere and the other in a vacuum 
chamber. Obviously in these two experiments the 
values of F are very different, i.e., F(patm) 7^ Fq 
where Patm is the energy density of the earth at- 
mosphere and so one may expect that such an ex- 
periment could distinguish between Palatini theory 
and general relativity. However, according to our 
above analysis, after coarse-grained averaging the 
relevant value of F in the atmosphere is nothing 
but Fq because NVp ^ Vtot- This means that the 
above thinking experiments will not work. 

So, in conclusion, at classical level the motion of par- 
ticles, cosmology and astrophysics in Palatini modified 
gravity theories are indistinguishable from the results of 
general relativity plus a cosmological constant. However, 
it must be emphasized that these are not equivalent the- 
ories. As we mentioned above, the internal structure of 
a particle in Palatini theories is generally different from 
that in general relativity (though we cannot measure the 
differences in masses classically) . This point will become 
relevant when we consider the structures of atoms and 
neutron stars in these two types of gravitational theo- 
ries (see § IVI|) . Furthermore, in Palatini theories the 
evolution of a region of space where the natural micro- 
scopic energy momentum tensor, T°f,, truly was that of a 
continuously-distributed pressureless dust would be dif- 
ferent from that of a region with the same coarse-grained 
density but where the matter was microscopically con- 
tained in small particles. 



V. THE ELECTROMAGNETIC FIELD 

In the previous section we have considered the dynam- 
ics of classical particles in Palatini theories. In this sec- 
tion we investigate how radiation, especially that due to 
the electromagnetic (EM) field {i.e., photons), behaves in 
these theories. This is particularly important since it is 
necessary to understand the propagation of light in order 
to correctly interpret cosmological observations. 

We start again from the microscopic theory and will 
be careful about the averaging procedure. While in the 
previous sections we have modelled the classical particles 
as occupying tiny portions of the space and in between 
of them there is vacuum, the electromagnetic field per- 
meates in the space diffusively and can be treated as 
a continuum. Consequently its averaging procedure is 
somehow like the naive one and slightly different from 
what we have met for classical particles. 

In the (jab frame the energy momentum tensor of the 
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EM field is: 

where Fab — 2d[aAi,]. We begin by considering the simple 
case of f{TZ) theories. In these theories, with ^ — TZ, 
Eq. ([6]) gives: 

2/($) = (58) 

so <J> takes its constant vacuum value and does not depend 
on Fed, and Eq. ([3]) gives: 

and so g"^ = Z.^.g"''- We also have that = j2^. Thus 
the effective energy momentum tensor, T°6, that appears 
in the gat-frame microscopic Einstein equation, Eq. ([7]), 
is: 

Kem = -g^'^g'^'F^dFeb (59) 
+ ls\g^fg''^FedFef - ^-fb- 

and F($)//^ — const. Apart from the addition of a 
vacuum energy term {—V{^)/f%{^)6^), the effective en- 
ergy momentum tensor in the gab frame has precisely the 
same form as it does in the gab frame. It follows that in 
vacuum: 

VaF"'' = 

where ^ agbc — and we have raised the indices of F°-^ 
using gab- The vacuum energy term generates an effec- 
tive cosmological constant, but other than that, the com- 
bined system the Palatini microscopic vacuum Einstein 
Maxwell equations in both metric frames is precisely the 
same as it is in GR. There is therefore no need to con- 
sider further the effect of averaging in these theories, or 
the way light propagates, in f{TVj theories since the they 
will be precisely as they are in unmodified GR. This 
comes about because in these theories the two metrics 
are conformally related, and the electromagnetic action 
is conformally invariant. As such the energy-momentum 
of the EM field does not source the extra Palatini degree 
of freedom encoded by <&, and the evolution of the EM 
field is not affected by the modification of gravity. 

More generally, however, in fiJZ, W'^Rab) the two met- 
rics gab and gab would instead be disformally related. In 
these theories then the microscopic behaviour of the elec- 
tromagnetic field will be altered, and hence its large scale 
behaviour may also change. It is most straightforward to 
see how the EM field behaves on macroscopic scales if we 
would in a frame where gravity is certainly no stronger 
than it is in GR over microscopic scales i.e. the gab frame. 
We therefore consider the form of the energy momentum 
tensor, T"f,, in the gab frame. Because of the complicated 
relationship between and g"'', as given by Eqs. ([31 
it is not possible, in general, to write down a simple 



expression, in terms of Fab, for T°f,. We therefore con- 
sider the relatively simple, but cosmologically interesting 
case, where the EM field is microscopically disordered 
and hence describes a bath of radiation. 

Choosing an appropriate coordinate system [t,x,y,z) 
so that at some point (to,xo), gab — Vab, we may write 
the components of this T"-^ explicitly, for example, 

= EyHz — HyEz 
i=y,z 

T^'^ = —{ExEy + HxHy). 

These are respectively the energy density, heat fiux, pres- 
sure and anisotropic stress terms, in which E, H are the 
strengths of the electric and magnetic fields. 

Since a totally disordered electromagnetic field neces- 
sarily, on average, has no preferred direction, when mi- 
croscopic fluctuations in the field are averaged over we 
must have: 

{El) = {ED {El), {Hi) = {Hi) = {Hi) 

Also 

{EyH, - HyE,) = [y^x] = {z^x)= Q 

and 

{E^Ey) = {EyE,) = {E^E,) = 0, 
{HxHy) = {HyHz) = {HxHz) = 

because of the lack of phase relations amongst the differ- 
ent components of the field strengths in the disordered 
EM field. The above results imply that on macroscopic 
scales T°'b for a disordered EM field behaves like a perfect 
fiuid with p = 3p and vanishing heat flux & anisotropic 
stress. This is precisely how EM fleld is treated cosmo- 
logically in standard GR. However, in Palatini theories, 
only represents the energy momentum tensor of the 
(jab frame. There gravity can be very strong over small 
scales, and hence behaves in a significantly different fash- 
ion to how it behaves in GR. We prefer to consider the 
gab frame energy momentum tensor, Tab, since in this 
frame gravity behaves no differently to how it behaves in 
GR. 

How the averaging actually works in the case of Pala- 
tini fiJZ,TZ'^^Rab) gravity depends on the microscopic 
field equations Eq. d?]), and here comes the difference 
between the fiJV) and f{TZ°'''Rab) cases. For convenience 
we decompose the symmetric tensors Rab and Tab respec- 
tively as 

Rab = Kab = ^UaUb + Sfafc + 2-U(aTb) + T,ab, (60) 
Tab = pUaUb + P^ab + 2Ui^aqb) + '^ab, (61) 
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where Ua is the four velocity of the observer {g^'^UaUf, = 
Cab = gab + UaMfe is the projection tensor to the 
hypersurface perpendicular to m, and p,p, qa, TTab are the 
energy density, isotropic stress, heat flux and anisotropic 
stress. Note that u^^ab — u°'qa — u°"'^ab — and TZ = 
— A + 3S. Throughout we use the convention that = 
g^^R^b but iC^b = r'K,b = r^Rcb- 
We note that Eq. ^ gives: 



R - Rabg"" - ^ {f.^r' + V.xK"') , (62) 



(63) 



Q2 Q2 

where the last two equalities follow from Eq. ^ and the 
definitions: $ = Kabg"'', X = KabK"''. Eq. g]) then gives 
(where all indices are raised w.r.t. gab)- 



up 



1 



A(/,$ - 2/,^A) + 2/,xT,T- + -/ 



(64) 



^«ga = (/,$ + 2/,x(S- A))Ta + 2/,;,T^I],,, (65) 
= (/,$ + 2/,xS)S - ^ + ^ (SafcS'^^ - T.Te^q) 

(67) 



-— :^4afc (^^cd^ — lei 



We also have: 

$ = as -A, (68) 
X = A2 + 3S2-2TcT^ + I]"6S''a, (69) 

and from Eq. ©: 

/,*$ + 2/.xX-2/ = K(3p-p). (70) 

We are concerned primarily with the form of the effec- 
tive energy momentum tensor Tab which is defined by 
Eq. (O. Defining Ua = \ua and ^ab ^ gab + UaUb so that 
UaUbg'^'' 



-1 i.e.. 



g2 _ 1 

we decompose Tab thus: 

Tab ^ pUaUb + pCab + 2U(^aqb) + T^ab (71) 

where (raising indices with g"'') we have If^qa — 0, 
U^^ab = and g^'^^ab = 0. 



We note that 



1 



«;T% = — (r%~F($,x)<5%) 



and so: 
P 

P 



_P_ _ 2a2/^ 

P 2a2/. 
2- 
A 



3Q- 



''TV- 



/,^<1> + 2/.^x - / 
2AtQ2 

(/,^$ + 2/,^x - ./) 
2«;Q2 



9" = 02 + 2/,x2g" + 2/,xS%g^] (72) 



Q2 



A similar expression to the above can be found for Sah 
although for our purposes it shall be not be needed. We 
note that if the average values of qa and Eab vanish, then 
so do the average values of qa and T^ab- 

The symmetries of (T"{,), i.e., {qa) = (TTafc) = imply 
that {qa) = {TTab) = 0, and that: 



(k/5 — 3kp) — 



2il^<f + 2/,^x - /) 



and by Eq. ^ wc have /,$$ + 2/,^x = 2/ so: 



(73) 



The average form of f°'b would only be equivalent to ra- 
diation with some cosmological constant if its trace (as 
given by the above equation) was constant, ft is evident 
from the above equations, that even if we (incorrectly) 
replace any quantity Q with {Q) in the above equations, 
e.g., TTab {TTab), ^a'^b — > (Ta)(Tb).) and usc the av- 
eraged values {p) = 3{p) and {qa) = {TTab) — 0, then we 
would find {Ta) = {^ab) = 0, and conclude that all av- 
eraged quantities involving or "Eab also vanish. If we 
accept this as true, we still have by Eqs. ([64]) and ((66)) : 

1 



{A)=H^{{^p--f)), 



i^)-H,i- 



up 



\f)) 



(74) 
(75) 



where 



H±{x) 



4/,, 



^ f% 
16/2 



± 



2/, 



It is then clear that generally when /^^ ^ 0, ($) = 
(3S - A) 7^ const. Similarly (x) = (A) ^3 (S)^ ^ const, 
for p 7^ 0. Hence unless = 0, we have that (/($,x)) 
and then also (i?) depend on p and hence on p. This 
means that in general f{TZ,TZ'^''Rab) theories (in contrast 
to the f{TV) models) a radiation dominated Universe will 
not obey the expansion law a (xt^, and moreover it will 
not even behave as a Universe with radiation and some 
cosmological constant. These deviations will be particu- 
larly significant for early time cosmology if f{TZ, TZ'^^Rab) 
is chosen to make the deviation from GR more significant 
at high densities. This represents the first complication 
associated with the behaviour of the EM field in Palatini 
theories. 

In deriving the above, we assumed incorrectly that be- 
cause (Ta) = (Sab) = 0, we could take the average values 
of quantities such as TaT° to vanish as well. Even though 
{qa) = 0, in general for a disorder electromagnetic field: 



{K^-qa) ^ 0{E^H^)k' 



and 



0{E'' + H^)k^ 
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are obviously nonzero. It then follows from Eqs. 
Wl\ that Kp, Kp, n'^q'^Qa, K^n°'''TTab and K^q'^q^iTab are re- 
lated to the five (scalar) unknowns A,S, T^Tq, Saf,E"'' 
and T^T^Safc. Also unless qa = T^ab = 0, we do not gen- 
erally have TqT" = 0, E^bE"'' or TaT^T^ab- Indeed since 
TaT°, EabS"'' > 0, the average values of these quantities 
could be non- vanishing even though the average values of 
qa, T^ab, and Eafc do vanish. 

Now in application to cosmology suppose we have / 
7^ A ~ S (~ means of order, / ~ 7?, simply means 
that the correction term is significant), FA ^ FS ^ 1. 
Then it could be seen that (A) ~ (5) ~ 0{kp) and 
(T.T'^) ~ - 0(«:V'), so that 



(A^ 



/■^2\ 



and the terms involving TaT", Y.^'^Yjab m Eqs. ([M1 - IS7)) 
are equally important as other modification terms. So 
when calculating the cosmology in a radiation dominated 
Universe we must take them into account. This is the 
second complication of the Palatini f{TZ'^''Rab) model. 
We would still, however, find that generally (i?) depends 
on p, signalling a deviation from the GR behaviour, when 
Ix ^ 0. 

These two complications make the Palatini f{TZ'^''Rab) 
model less trivial than its f{Tl) counterpart: for the latter 
the cosmologies for both radiation and matter dominated 
Universes are the same as in ACDM, while for the former, 
as we referred to in above, the radiation dominated Uni- 
verse could behave rather differently from ACDM, even 
violating the a oc 1 2 law of expansion. It will be interest- 
ing to consider such models in more details, but due to 
its complexity this is beyond the scope of this paper and 
will be further investigated elsewhere. 

But how does this difference from Palatini f{TZ) grav- 
ity and from the behaviour of classical particles arise? In 
the case of classical particles, there are no interactions 
between separated particles other than gravity and this 
is why the cosmology is like that of GR with a cosmolog- 
ical constant; the modification to GR simply alters the 
internal structure of the particles. The radiation field, as 
we said above, could be treated as a continuum and the 
new 'interaction' due to the modification to GR (which 
acts at a point!) exists everywhere; so photons feel the 
modification everywhere and the corresponding cosmol- 
ogy is changed: this is essentially the argument of a naive 
averaging. The Palatini f{Tl) is immune to this effect 
because in this model the interaction depends on p — 3p 
which is zero identically. In general f{TZ, TZ'^^Rab) models 
the interaction does affect the propagation of photons. 



VI. CONSTRAINTS FROM ATOMIC PHYSICS 

In this §we provide a particular example of how micro- 
scopic physics tightly constrains the properties of Pala- 
tini theories (see also [i^l). We consider how the Palatini 
modification alters energies of the photons that are emit- 
ted when an electron transitions from one atomic energy 



level to a lower one. This analysis is particular to f{TV) 
theories. A similar calculation could be performed for a 
general f{TZ,TZ'^^Rab) theory, but this would be signifi- 
cantly more complicated due to the fact that the light 
propagation is also altered in these theories. 

Consider the total action for a charged fermion field 
vpp: 



Stat — Se 



rav ~ 



where 



Srn — 




(76) 



(77) 



gA'^x {^Y (e^Va - mo - iqg^Aa) *f 



-ab-cd 771 771 
g g faci'bd 



(78) 



where Aa is the electromagnetic field, q is the charge of 
the particle, f?" are the curved space-time analogue of the 
Dirac 7 matrices and q°'q^ -f q^q'^ = 2^"''. 

We will perform our calculations in the Einstein frame, 
i.e., where the metric is gab- We do this because, as we 
noted above, in this frame the local curvature of space- 
time is certainly of a similar magnitude to that which 
one would expect from general relativity. As such the 
approximation gab ~ rjab in the frame of a laboratory ex- 
periment will be equally as valid in Palatini theories as it 
is in standard general relativity. The same is not neces- 
sarily true for the Jordan frame where gab is the metric. 
Working in the Einstein frame should be viewed merely 
as a computational convenience and should not be viewed 
as attaching any special physical meaning to this frame. 
The observable quantities we will extract from our cal- 
culation will be independent of the choice of frame. 

We now convert the action, as given by Eq. (|76|) . to 
the Einstein frame. In this frame we have gab — f,igab, 
= R, ^ 



f%V-9, and so 



-gf{<i>) = V^/*<i> + V^ [/(*)-/,*$] 

^ ^[R-2kV{<S>)], (79) 



where again ^ = TZ and V{^) = (/,*$ - /) /2k. Thus 
Eq. ((75|) becomes: 



•Stot — 



C I c(eff) 



where 



c 

^gra 

'^m '-'77 




JV-gd.j 



(80) 

(81) 
(82) 



The total energy momentum tensor associated with vari- 
ations in Sf-^f with respect to gab is conserved. 
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We rewrite S'm in the Einstein frame as follows: we 
define new fields 



and then: 



and 



so that 



-1/2 a 



We then have: 



(83) 



in which Va satisfies V aQbc = and m($) = mof ^ ■ 
Varying this action with respect to tp, we arrive at the 
modified Dirac equation obeyed by electrons in these the- 
ories: 



Varying the action with respect to $ gives us: 
2/($) - 



f2 



Km{^)(fip. 



(84) 



(85) 



With the equations written in this form, the effect of 
the Palatini modifications is manifest: it is to make the 
mass of the electron, m(<i>), a ^-dependant quantity (see 
e.g. [13, m, |4§| for considerations of other local density 
dependent quantities). Since the mass of the electron 
depends on $, it also depends on the local density of 
matter, and hence also on the local electron density. In an 
atom the peak local electron density generally decreases 
as the energy of the orbit decreases, and so the effective 
mass of the electron will be different for different energy 
orbits. This will lead to potentially detectable deviations 
from the standard model of particle physics. We now 
quantify these deviations. 

We wish to consider the energy eigenstates of an elec- 
tron orbiting an a hydrogen nucleus. We therefore take 
if oc exp(— ii^ei) where is the electron energy, and for 
the hydrogen atom qAq = a-Eu/r = — V{r) where aEM is 
the fine structure constant. From Eq. the electron 

obeys: 

Ac(y5 = -lakOkf (f + m{9)fiip, 

r 

where ai/3 = —f3ai, = 1 and a^jaj-) = Sij. We define: 



f F{x) 
^ = ( ^G{x) 



(86) 



- m($) 



ttEM 



Ec + m($) 



OEM 



F = cr • VG, 

G=-a-VF,, 



(87) 
(88) 



where Ui are the Pauli matrices. E is the total energy 
of the electron and includes the contribution from the 
electron rest mass, we therefore write E^ = toq + E'^. 

We calculate the modifications to the energies of the 
photon that is emitted when an electron changes from 
one energy level to another by assuming that the Palatini 
modification is small enough so that we may write /(<&) — 
<i>[H-e($)] where e(<I>) <C 1. Not aU Palatini theories can 
be written in this may, but considering we will find that 
any changes in e($) are constrained to be very small, we 
believe that is it highly unlikely that any theory, which 
cannot be written thus, could be experimentally viable. 

Assuming this form for /($) we have 

to($) = mo/v/Z*« mo [1- ^ (£ + $£,<!.)] EE mo+(5TO($), 

where (5TO($)/mo ^ 1. We solve Eqs. (|571[551) perturba- 
tively in the small parameter (5m($)/mo. 
To 0(((5m/TOo)") we have 

F — F, G = G and E'^ = Ec — a|]y[r7iof (aEM) 
where having defined y — x / a^, — 1/aEMmo: 



ttEM 



O I ^2 c , "EM 

\y\ 



F = ff-V(^)G, (89) 
G = -aEMCT • V(j,)F, (90) 



thus G ^ 0{aEMF). The unmodified energy levels, Ec 
are just what they would be in the standard model with 
an electron mass mo . We do not repeat a full calculation 
of those energy levels here. For the purposes of calcu- 
lating the large contribution perturbation to the energy 
levels we need only work to leading order in c^em • To the 
leading order in oemj F satisfies the Schrodinger equa- 
tion: 



F{l + 0{al^)). 



(91) 



To this order, for a state with energy level n and angular 
momentum {Im) we have to leading order in a-^M'- 



£^£n = -1/2^2(1 + 0(a|j^)) 



and we write 



where Yim{d^ 0) are spherical harmonics, r is the distance 
from the nucleus and 9 and cj) are angular coordinates. 
We normalize Fnim so that: 



d-'xF^i^ = 1. 
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To leading order in both oem and Sm/rriQ the electron 
number density, rie, is given by rig — and, so to this 
order by Eq. ((85|) . where for a state with energy level n 
and angular momentum (Im), we have $ = and 



To O^Sm/mo) we write 

F = F + SF, G = G + SG and 
We therefore have: 



(92) 



+ Am. 



(5™ G 



mo QfEM 



Am + (5m 








" m) 


mo 







-2a 



EM 



Am — Sm 








F + 




mo 







(93) 



and so SG ^ 0{aEMSF). To leading order in q;em then 
we have 6F — XF / a^j^ where: 



V,(|FpV(,)A) = 2 



Am — Sr 



mo 



(94) 



integrating this expression gives, for a state with an en- 
ergy level n, and angular momentum (/m): 



d^X(5m„i™|F„;™| 



mo 
2 



(95) 

where Sninim — Sm{^nim) and we have defined: 



iQ)nln 



d^x(5(a;)ne („;m)(x). 



In a local inertial fi'ame we ignore the energy stored in 
the gravitational field, however we can see from Eq. ([5^ 
that even when gab — rjab, the effective matter action in- 
cludes a contribution from the effective potential V{^). 
Only when the contribution from V{^) is included is the 
energy and momentum conserved with respect to gab- 
Thus (excluding the energy of the nucleus which we as- 
sume to be independent of the electron energy level), the 
total conserved energy in this set-up is given by: 



E-y 



where E-y is the energy stored in the photon field and 
when gab = rjab- 



E^ 



= d 



^3. Vi^) 



(96) 



Assuming as we have that /($) = $(1 + e($)), we have: 



J $ 



2k 



Thus taking $ = $nim we have: 
mo 



E 



v(<s>) 



Im 



(97) 



where we have used Eq. ([92)1 . Thus: 

E^ + Ev(^) ^mo + E^{nlm) - ^ (e)„;„ (1 + 0(q|m)). 

Since Etot is a conserved quantity in this frame, when an 
electron changes from a level with quantum numbers (nl) 
to a lower energy level with quantum numbers (ji'V), a 
photon with energy, E^ must be emitted where: 



E^ =^e(nl)~^e(n'l') 
where we have defined: 



(98) 



1 ' 



m— — l 



We define _E^(nl, n'l') be the energy of the photon to ze- 
roth order in e. We then liave: 



1 + A"?' 



^7(nl,n'l'), 



(99) 



where, using _Ec(nl) = —a'^^mo/2n'^ to leading order in 
ttEM, we have to leading order in aEM: 



A n'l' _ 
^nl — 



/n'l' 



(n 2 — ri' 2) 



(100) 



'EM 



If A were the same for all transitions then this could be 
account for simply by a slight alternation of the electron 
mass. To constrain A we consider therefore how the en- 
ergy of a photon released due to one transition changes 
relative to that emitted in another transition; this is in- 
dependent of the electron mass i.e. we consider the ratio 
of E^ for one transition with that for another. It is also 
independent of one's choice of frame. If t is the time in 
a local inertial frame in the Einstein frame, and i the 
time in a local inertial frame (LIF) defined with respect 
to the Jordan metric, then, since gab = f,<s>gab, we have 
= /,*dt. Thus if a photon has energy E^ in LIF of 

1/2 

the Einstein frame metric, gab, it has energy E^f ^ in 
a LIF of the Jordan frame metric, gab- If the energies of 
two photons are measured at the same place, the ratio 
of those two energies is a frame independent quantity, as 

1/2 

the frame dependant scaling factors (/^ ) cancel. Thus, 
even though we have performed our calculation in the 
Einstein frame, the quantity which we will constrain is 
independent of this frame choice. 

We now consider two transitions for where measure- 
ments have been shown to agree to the standard the- 
oretical prediction to a high accuracy [i^]- Firstly we 
have the 151/2 ~ 2S'i/2 transition. For this transition 
(nl) = (20) and {n'l') = (10). Secondly, we consider the 
2>S'i/2 — 81)5/2 transition. For this transition (nl) = (83) 
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and {n'l') = (20). Thus: 

^-y(83;20)i^-^(20;10) 

64 „ , 



/83 



28 

5q:eM 



16 
21 



/20J 



/83 



/20 



'lOJ 



21 



(101) 



The measurements of Ref. 42| provide the foUowing con- 
straint: the magnitude of the left hand side of the above 
equation to be smaUer than 8 x 10^^° and so: 



16 , ^ 
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< 8 X 10"^'^. 



(102) 



This represents a very strong constraint on the properties 
of Palatini f{Ti) theories. In particular we can see that, 
writing /($) « $(1 + e($)) when ^/k is of the order 
of the electron cloud density in hydrogen, then changes 
in e($) are constrained to be very small. If we expand e 
about some appropriate value of $ and find that to linear 
order we have 

e($) « const + eo^/bH^, 

where Hq is the value of the cosniological constant today, 
Eq. (|102p gives the very strong constraint: 



\eo\^\f"{mo/f'm<^^w 
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VII. DISCUSSIONS AND CONCLUSIONS 

Much of our intuition about how the microscopic be- 
haviour of gravity affects physics on large scales is based 
upon Einstein's general relativity. In this article we show 
that such intuition cannot simply be applied to modi- 
fied gravity theories without a detailed analysis of the 
energy-momentum microstructure. Indeed, naively aver- 
aging over the microscopic structure will generally lead 
one to make incorrect predictions, and inaccurate conclu- 
sions as to the validity of the theory. In particular, the 
nai've averaging procedure is invalid in Palatini theories. 

For classical particles, we show that the relative mo- 
tion of particles in Palatini theories is indistinguishable 
from that predicted by GR plus a cosmological constant. 
This means that the cosmology and astrophysics (except 
in some extreme environments such as neutron stars) of 
Palatini f{TZ, TZ°^^Rab) models are identical to that of GR 
plus a cosmological constant. The above result can also 
been shown using our correct averaging procedure, as 
given in ^ilVEl This particularly means that the fine 
tuning problems associated with the cosmological con- 
stant are not alleviated in Palatini theories. It should be 
stressed that although the Palatini theories predict the 
same cosmology and astrophysics as GR, they are com- 
pletely different theories: not only because they predict 



different internal structures of particles, but also because 
they behave differently from GR in the presence of elec- 
tromagnetic field and in the atomic physics. 

When coming to electromagnetic fields in Palatini the- 
ories, things becomes a bit different. In contrast to clas- 
sical particles, which are tiny clumps of energy density 
in between of which there is vacuum, the EM field per- 
meates in the space and the nai've averaging actually 
works. However, when performing the averaging one 
should also take into account the fact that the field equa- 
tions are microscopic and that at microscopic level the 
EM field is random and disordered. This make the Pala- 
tini f{TZ°'^Rab) model less trivial than its f{Tl) counter- 
part: for the latter the cosmologies for both radiation and 
matter dominated Universes are the same as in ACDM, 
while for the former the radiation dominated Universe 
could behave rather differently from ACDM, even violat- 
ing the a (X t2 law of expansion. This difference from 
Palatini f{TZ) gravity and from the behaviour of classi- 
cal particles arises because, in the case of classical parti- 
cles, there are no interactions between separated particles 
other than gravity. This is why the cosmology is like that 
of GR with a cosmological constant. The modification to 
GR simply alters the internal structure of the particles. 
The radiation field can be treated as a continuum and the 
new 'interaction' due to the modification to GR (which 
act at point!) exists everywhere. So photons feel the 
modification everywhere and the corresponding cosmol- 
ogy is changed. The Palatini f{Tl) is immune to this 
effect because in this model the interaction depends on 
p~3p which is zero identically. In general f{TZ,TZ"'''Rab) 
models the interaction does affect the propagation of pho- 
tons. In summary, the Palatini modifications can affect 
the propagation of photons (EM field) and even change 
the cosmic expansion during radiation domination. Ob- 
servational data on Big Bang Nucleosynthesis could then 
place some constraints on these models. 

Interestingly, although Palatini f{TZ) theories were de- 
signed to modify gravity on large scales, they actually 
modify physics on the smallest scales {e.g., the energy 
levels of electrons) while leaving the larger scales practi- 
cally unaltered. We show that the observables in atomic 
physics, such as the energy levels, can be very sensitive to 
the Palatini modification to GR, and indeed experimen- 
tal data places extremely stringent constraints on any 
deviation from GR. In general, before considering any 
astrophysical consequences of a modified gravity theory, 
it is important then to check that it does not make un- 
realistic predictions for atomic physics. 

One may wonder why the same averaging problem does 
not appear in the metric f{R) gravity models. The an- 
swer is that, within the metric approach, averaging over 
microscopic scales is generally no less straightforward 
than it is in GR; this is because in both cases all degrees 
of freedom are dynamical. These dynamics generally en- 
sure that the field equations, for all degrees of freedom, 
are approximately linear for small-scale structures. In 
contrast, averaging in Palatini models is not so trivial 
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since the new degree of freedom is non-dynamical, and so 
its field equation remains non-linear even on the smallest 
scales. This non-linearity introduces an averaging prob- 
lem that is specific to Palatini theories [52*1 . And this 
explains why the cosmological behaviour of these theo- 
ries can be very different from what has been suggested 
in the literature [H, O, Sll . 

It is usually the case that modified gravity theories, 
other than those considered here, predict extra non-linear 
terms in the microscopic Einstein equations. It is there- 
fore both important and interesting to check that av- 
eraging, both over microscopic and cosmological scale 
structures, does not significantly alter the macroscopic 
behaviour of the theory from that which might have oth- 
erwise been expected from the microscopic equations. As 



much as the importance of any back-reaction from aver- 
aging remains an open problem in General Relativity, as 
we have illustrated in this work, it is likely to be an even 
more significant problem in modified theories of gravity. 
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